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Abstract 

We investigate the structure of the invariant measure of space-time chaos by adopting an "open- 
system" point of view. We consider large but finite windows of formally infinite one-dimensional 
lattices and quantify the effect of the interaction with the outer region by mapping the problem on 
the dynamical characterization of localized perturbations. This latter task is performed by suitably 
generalizing the concept of Lyapunov spectrum to cope with perturbations that propagate outside 
the region under investigation. As a result, we are able to introduce a "volume" -propagation 
velocity, i.e. the velocity with which ensembles of localized perturbations tend to fill volumes in 
the neighbouring regions. 
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I. INTRODUCTION 



Since the discovery of deterministic chaos, it has become clear that unpredictable be- 
haviour can not only be the outcome of a stochastic dynamics, but also of a few nonlinearly- 
coupled degrees of freedom. Even more important is to notice that these two classes of 
behaviour can be distinguished without making any assumption on the underlying model. 
This is made possible by embedding a supposedly recorded time series Uj = u(t = jr) 
(where r is the sampling time) into a space of dimension L (i.e., introducing the vec- 
tor \jf ] = {uj ,Uj + i, . . .Uj+L-i}) and thereby estimating the fractal dimension D(L). In 
stochastic processes, D(L) = L, since the variables at different times are mutually indepen- 
dent (at least below a certain threshold, which clearly depends on r). In low-dimensional 
chaos, when L is increased, D(L) saturates to a finite value because of the functional depen- 
dence among the variables. This is the starting point of nonlinear time-series analysis, an 
ensemble of tools, developed in the last years to reconstruct a deterministic model starting 
from raw data,Q]. 

More subtle is the difference between stochastic and deterministic signals, when the latter 
ones arise from a high- dimensional dynamics (as in space-time chaos), since D(L) increases 
with L in both cases. In this context, it is still possible to distinguish between the two classes 
of behaviour, provided that a more refined analysis, based on the effective, or coarse-grained, 
dimension D c (e, L) is developed, where e represents the resolution of the coarse graining. 

In fact, in stochastic systems, all variables turn out to be mutually independent below 
a fixed threshold e s that is independent of the embedding dimension L, so that D c (e, L) is 
essentially equal to L for e < e s . On the other hand, in space-time chaos, it has been conjec- 
tured that new degrees of freedom appear only at the expense of progressively decreasing the 
observational threshold. A loose explanation for this behaviour is based on the observation 
that the degrees of freedom associated to far-away regions are almost decoupled from the 
evolution in the observation point |2j. Although this statement looks very reasonable, it is 
not at all easy to substantiate it with solid arguments. Indeed, very little progress has been 
made in the last decade and what is known is still mostly based on heuristic arguments. 

In this paper, we attack the problem of characterizing space-time chaos by constructing 
a spatial rather than a temporal embedding, i.e., by referring to the hypothetical signal 
Uj = u J (t), where the time t is large enough to ensure convergence to the attractor, while j 
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now labels the lattice sites (for the sake of simplicity we limit ourselves to considering one- 
dimensional lattices). Accordingly, the effective dimension D c (e,L) now counts the number 
of degrees of freedom that can be resolved, with resolution e, in a window of length L 
embedded in a supposedly infinite system. This problem is conceptually equivalent to the 
previous one, the main difference being that time and space axes have been exchanged. It is 
precisely this difference that allows us using the more standard tools developed to investigate 
the invariant measure of chaotic systems. 

In the investigation of space-time chaos, closed systems of finite length L are typically 
considered. In such a context, the classical tools for the characterization of chaotic dynamics 
such as the Lyapunov exponents describing the evolution of infinitesimal perturbations, can 
be effectively implemented. As a result, it has since long been recognized that a limit 
Lyapunov spectrum does exist for L — > oo, thereby inferring (from Kaplan- Yorke and Pesin 
formulae) the extensivity of space-time chaos |3j: in fact, both the fractal dimension and 
Kolmogorov- Sinai entropy are proportional to the system volume (length, in one dimension). 

At variance with this approach, here we adopt an "open" -system point of view, i.e. 
the finite window of length L is part of an infinite system whose evolution is taken into 
account as well. In a sense, the two methods are reminiscent of the microcanonical and 
canonical ensembles of statistical mechanics; the open-system approach is indeed a possible 
way (though not the most effective one) to perform canonical simulations of Hamiltonian 
systems. Unfortunately, at variance with equilibrium statistical mechanics, here there is 
not a prescription such as the Boltzmann weight to estimate, a priori, the probability of 
each configuration. All we have at our disposal is the nonlinear dynamical law with the 
additional difficulty (in comparison to the "closed" -system approach) of having to deal with 
perturbations propagating from the outer to the inner region and viceversa. 

Collet and Eckmann have developed a method that more than any other has inspired 
us in developing the approach outlined in this paper. Their idea is that the structure of 
the invariant measure in a window of fixed length L at a given observational scale can be 
described by following in time the convergence of a suitable ensemble of initial conditions 
towards the attractor. In fact, such an evolving "cloud" of points provides a natural covering 
of the attractor that becomes increasingly sharp upon time evolution. Turning the time 
dependence of infinitesimal ellipsoids into a dependence on the observational resolution e is 
a possible way to explain the Kaplan- Yorke formula in standard finite dimensional systems. 
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However, the extension of such ideas to open systems requires that the very concept of 
Lyapunov spectrum be revisited so as to describe perturbations that evolve also outside 
the region that is currently monitored. We shall see that a meaningful characterization of 
the perturbation dynamics can be obtained only by realizing that the problem involves two 
scaling parameters that must be simultaneously let diverge to infinity: the window length L 
and the evolution time T. In fact, we end up introducing a Lyapunov spectrum that, besides 
depending, as usual, on the ratio of the label of each exponent by the length L, depends 
also on T/L. This is not simply an additional technical difficulty, but the key element that 
allows mapping the characterization of pertubation evolution onto the characterization of 
the invariant measure in spatially extended systems. In fact, expressing the variable T in 
terms of the observational resolution allows sheding further light on the dependence of the 
effective dimension D c (e, L) on both L and e. While our analysis does confirm the functional 

n n 

dependence conjectured in |5|, |6j, we find different expressions for the coefficients. 

A brief summary of the known results is presented in Sec. II. Sec. Ill is devoted to 
a detailed justification of the Kaplan- Yorke formula by following an approach that can 
be most easily extended to open systems. A relevant result of our analysis is the above 
mentioned extension of the concept of Lyapunov spectrum: this issue is the core of Sec. IV, 
where we also illustrate the implementation of the method in various classes of coupled map 
lattices. Finally, in Sec. V, we discuss the implications of the Lyapunov analysis on the 
fractal properties of the invariant measures in open systems and briefly discuss the further 
corrections expected to arise from the boundaries. 



II. THE STATE OF THE ART 



As already anticipated in the introduction, we aim at characterizing the structure of the 
invariant measure of space-time chaos. For the sake of simplicity, most of the analysis will 
be restricted to (one dimensional) lattice systems. A detailed description of the scaling 
properties of a given set is contained in the "effective" dimension 

„ , . dH , _ . 
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where H{e) is the entropy of the set covered with boxes of size e, H = —J^Pi hipi, p% being 
the probability of each box 1 . As shown by Renyi long ago, different definitions of entropy 
can be given, by replacing the logarithmic average of Pi in the H expression with averages 
of its moments. As, in general, there are tiny differences among the various entropies, we 
will always refer to H without specifying which average is being taken. 

To our knowledge, the concept of a resolution-dependent entropy was first introduced 
by Kolmogorov and Tikhomirov 3], who called it ^-entropy In principle, an e-dependent 
dimension is an ill-defined concept, since it is not invariant upon change of space parametriza- 
tion. For this reason, one has to take the limit e — ► 0. In fact, only in this limit, T> c (0) 
becomes a dynamical invariant that is strictly related to other dynamical invariants such 
as the Lyapunov exponents. However, besides the asympotic value T> c (0), also the possibly 
slow dependence on e at high resolutions can be universal, thus conveying meaningful infor- 
mation about the structure of the set of interest. This is precisely one of the reasons why 
^-entropies have been introduced to quantify the cardinality of different classes of functions 
(such as, e.g., the entire functions) P,Q|. 

In the study of spatially extended systems, one is faced with the difficulty of accounting 
for the dependence on a further parameter besides e, namely the system size L. However, 
previous studies of closed systems have clearly revealed that, for sufficiently large L, the 
coarse grained dimension is still an extensive quantity: [3] 

D c (0, L) oc dL (2) 

where d can be interpreted as the dimension density (i.e. the contribution to the dimension 
per lattice site) that can be determined from the Kaplan- Yorke formula (see next section). 

In open systems, it is instead clear that D c (0, L) = L, since the infinitely many degrees 
of freedom ruling the outer part of the chain act as a sort of stochastic source ja in 
order to understand how the two results can be reconciled, it is necessary to investigate, in 
the case of open systems, the simultaneous dependence on both e and L (within the closed 
system approach, the resolution does not play an important role - see next section). As 
long as e and L are respectively small and large enough, the scaling dependence on the two 
parameters is expected to be universal. 



1 Rigorously speaking, one should refer to the optimal covering. We implicitly assume to have made such 
a choice 
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The first consistent conjecture about this problem was formulated by Korzinov and Ra- 
binovich 

vd 2 

D c (e,L) = dL \ne-A, (3) 

? n 

where a is again the dimension density, 77 is the Kolmogorov- Sinai entropy density [3J, v is 
the propagation velocity of disturbances, and A is a non-better-specified parameter. This 
equation parallels the analogous expression proposed by Tsirimg for the symmetric problem 
of the dimension of a scalar time series recorded at a single point j^. One can notice that 
Eq. © allows reconciling the apparently contradictory expectations for closed and open 
systems. In fact, if the limit e — > is taken before the limit L — > 00, D c /L diverges (though, 
in reality, it could not become larger than 1 - this inconsistency is due to the perturbative 
character of the above formula); if the order of the limits is reversed, then, D c /L converges 
to the expected (closed system) value d. 

Unfortunately, the derivation of the above formula depends on several assumptions that 
cannot be directly checked. Moreover, it is rather unlikely that more accurate numerical 
simulations will provide clean enough data to draw definite conclusions. It is therefore 
compelling to make some progress on the theoretical side, even at the expense of introducing 
strong simplifications. This is the route already undertaken in Ref. Q], where the limit case 
of weakly coupled maps has been considered. 

If the invariant measure is assumed to cover a linear subspace, it is possible to obtain a 
detailed description of it h|. In fact, under this approximation, it is possible to implement 
global methods such as singular-value decomposition (SVD) technique. In general, the 
usefulness of SVD is limited by the presence of nonlinearities that induce bendings which, in 
turn, do not permit extracting information about the local thickness of a given set. Let us, 
for instance, imagine a slightly bended segment embedded in a two-dimensional plane. SVD 
will tell us that the set of points belonging to the segment is characterized by two non-zero 
orthogonal widths even though the set is strictly one- dimensional. If one can restrict the 
discussion to linear subspaces, this problem does not arise and a global method like SVD can 
be effectively used to extract local information. In Ref. it has been assumed that the 
invariant measure of the infinitely extended system is the linear superposition of a subset of 
all possible modes in a given basis (e.g., all Fourier modes with wavenumber smaller than 
a prescribed threshold). The fraction d of the active modes is the dimension density of 
this given space of functions. The corresponding problem of characterizing the projection 



onto a space of finite width L can be addressed by computing the eigenvalues of a suitable 
correlation matrix. As a result, it has been found that the observational resolution e and 
the effective dimension D c are connected by the scaling relation 

\ne = -LF(D c /L), (4) 

where the function F{x) is identically zero for x < d, while it increases monotonously for 
x > d, starting with a finite slope. Solving the above equation with respect to x — D c /L 
and expanding the resulting inverse function, F _1 (hi£) for small values of hie, one finds the 
perturbative expression 

D c (e,L)=dL- l -^-^ E (\ne) 2 , (5) 

that has the same structure of Eq. @. Thus, Eq. (£Q) appears as its natural extension to 
arbitrarily small scales. 

The main drawback of this approach is the absence of any dynamics. In particular there 
is no way to link the function F(x) to dynamical invariants such as Lyapunov exponents. 

A more suitable starting point is represented by the study of the complex Ginzburg- 
Landau equation performed by Collet and Eckmann |4j. They have rigorously proved that 

D CE (e,L) = B L + B 1 /e 2 (6) 

represents an upper bound to D c . On the one hand, this formula confirms the existence of a 
leading term that is proportional to the system size. On the other hand, the above expression 
differs from the previous ones for what concerns the leading correction, that diverges faster 
than logarithmically for e — > 0. The question whether such a difference is to be attributed 
to the continuity of the space variable (and thus to the possibly larger number of degrees of 
freedom) or it is due to technical difficulties in improving the upper bound cannot be easily 
answered. On the basis of the results here presented we argue that, for lattice systems, the 
upper bound Eq. (0) can be improved, while for spatially continuous flows, the situation is 
yet unsettled. 

The idea behind the derivation of Eq. (JHJ) is the same that allows proving the Kaplan- 
Yorke formula for standard finite-dimensional attractors: given a set of boxes that cover 
the attractor, we can obtain finer coverings by simply letting each box evolve in time. We 
summarize the idea in the next section, as it will be useful for the generalization to open 
systems. 
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III. THE KAPLAN- YORKE FORMULA 



In this section, we discuss a method that allows extending the Kaplan- Yorke formula 
to open systems. It is both useful and necessary to start from the simple context of a 2d 
chaotic map (such as, e.g., the Henon map). Let us cover the attractor with a square So 
of size 0(1) (in general it will be a hypercube) and let us denote with St its image after 
t time steps. St provides a covering of the attractor at all times, even though stretching 
and folding transform it into a long and thin sausage. It is natural to divide St into boxes 
of size equal to its average width (for the sake of simplicity, we do not take into account 
multifractal fluctuations, that would not anyhow modify the following scaling arguments) 

e = 6 2 {t) = exp{A 2 t} (7) 

where A2 is the second, negative, Lyapunov exponent. This is the crux of our argument, as 
it suggests how time evolution spontaneously introduces an increasing resolution in phase- 
space. In fact, we can now estimate the fractal dimension T>(e) from the number of boxes 
N(e) needed to cover S t with resolution e, 

v(c Jj^n = _h^m = l _ A, 

In e In £ A 2 

where Ai is the positive Lyapunov exponent (in the philosophy of neglecting multifractal 
corrections, we do not distinguish between the positive Lyapunov exponent and the topo- 
logical entropy). S t provides a meaningful covering of the attractor only if the dynamics 
is invertible, otherwise the above equation would represent only a (possibly rough) upper 
bound. Eq. (|SJ) is nothing but the well known Kaplan- Yorke formula in 2d maps. It is worth 
recalling that although this derivation is rather sketchy, the equality is rigorous, provided 
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that T> is interpreted as the information dimension 

This approach can be extended to higher dimensional maps, but it requires an addi- 
tional assumption that seems to be generally valid, although not universally correct. Before 
discussing the most general case, let us first add a decoupled, contracting direction to the 
previous system, as this case helps clarifying the difficulties that arise in higher dimensions. 
The length of St remains unchanged, while its transversal section becomes an ellipse with 
two semi-axes of length 5 2 = exp(A 2 t) and <5 3 = exp(A 3 t), respectively (A3 being the addi- 
tional, negative, Lyapunov exponent). The question now consists in choosing the size e that 
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allows an optimal covering of St- In this special case, we know a priori that N(e) must not 
change, since we have not modified the attractor itself. If we choose e = 83, the resulting 
N(e) can be either much larger or smaller than before, depending on the relative size of 82 
and 83. The error of this choice is that the hidden structures not yet resolved at time t are 
contained only along the second direction and thus, we must fix e = 82- 

In this case, the right result has been obtained because, in a sense, we already knew 
the solution. Let us now refer to a generic finite-dimensional system and let N p denote the 
number of positive Lyapunov exponents. An initial hypercube S with edge- length 0(1) 
covering the attractor is stretched along the unstable directions and contracted along the 
stable ones. Since the attractor is bounded along all directions, the "excess" of length that 
is continuously produced along the unstable directions must be folded along the contracting 
ones. The crucial point is the assumption that folding generically proceeds from the least 
to the most contracting directions. Technically, this is equivalent to assuming that all 
contracting directions are filled until A& = Xw=i ^1 < (where the Lyapunov exponents are 
implicitly ordered from the largest to the most negative one), while the remaining most stable 
directions do not contribute (like the third direction in the above example). Accordingly, 
the "right" box-size to be adopted in the partitioning process of S t is 

£ = exp(A m t) (9) 

where m is the minimum fc-value such that A& < 0. From the corresponding number of 
boxes of size e needed to cover St, it is readily found that 

T> KY = m-l + £f=±. (10) 

This is the general form of the Kaplan- Yorke formula. 

In spatially extended systems of large length L, the Lyapunov exponent depends on the 
index / and L only through the scaling variable p = l/L, i.e. A; = A(p). Accordingly, by 
neglecting the fractional correction in Eq. (jlUj) . the dimension can be written as 

V KY = d KY L (11) 

where (Iky, implicitly defined by the constraint j^ KY dp\(p) = 0, can be interpreted as a 
density of dimension (see also Eq. Q) 0. 

All of the above discussion can be summarized stating that the number of boxes needed 
to cover an attractor with resolution e can be determined by letting a ball of size 0(1) evolve 
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until its width along the m-th direction (determined by the vanishing of A m ) is equal to e 
itself. In other words, Eq. Q is the core of the argument, as it allows transforming the 
dependence on t into a dependence on the resolution. 

In the above discussion we have implicitly assumed that each Lyapunov exponent is 
equal to its asymptotic value, independently of t (and thus of the resolution). As long as 
each Lyapunov exponent exhibit a slow dependence on time, all of the above discussion still 
applies, with the difference that the r.h.s. of Eq. (JTUJl depends on e (through the hidden 
dependence of the A's on t). In other words we see that the Kaplan- Yorke formula is ready 
to account not only for the asymptotic value of the fractal dimension but also for possible 
dependencies on the observational resolution. This is precisely what happens in the case of 
open systems. 

Before proceeding further in this direction, we need to introduce some notations: let xy 
and x_i_ denote two vectors defining the state variable on each lattice site, within, respectively, 
outside, the window of interest Wl- As in the previous discussion, the aim is to infer the 
fractal properties of the attractor from the density P(t, xy, Xj_) at time t (the initial condition 
P(0, xy, x_i_) being a constant distribution in a hypercube of radius order 1 that contains the 
attractor). The probability density can be usefully rewritten as 



where Q(t, xn, Xj_|xj_) denotes the probability density at time t conditioned to the initial 
state Xj_ of the external variables, while P±(x° ± ) represents their distribution. The integral 
over the "hidden" variables x^_ represents the first relevant difference with the previous 
case: the ignorance about their values contributes to dressing the probability density. We 
will discuss a bit this problem in the last section. 

Anyhow, even disregarding the effect of the integral in the above equation, the problem 
we have to deal with is more complex than the previous one. In fact, even if we consider 
initial conditions (xjj, x 1 }), (xjj, Xj_) that differ only inside Wl, the mutual difference does not 
remain confined to Wl, but rather spreads and propagate in the outer regions. Accordingly, 
we are faced with the problem of defining, in this context, the Lyapunov spectrum in a 
meaningful way. This is the goal of the next section. 

However, before discussing the generalization of Lyapunov spectra to open systems, we 




(12) 
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briefly present a heuristic explanation of Eq. (jUJ), since its derivation is based on the idea 
that a characterization of the invariant measure over increasingly fine scales can be obtained 
from the evolution of a set Sq covering the attractor. More precisely, in Ref. Wl has been 
split into two parts: the bulk B, where the effect of the external degrees of freedom can be 
neglected (over the time t), and the boundary b, where propagation must be, instead, taken 
into account. An upper bound to the number of boxes needed to cover the attractor has then 
been estimated as the product of the number of boxes needed to cover B, times the number 
of boxes needed to cover b. In the computation of both quantities one is faced with the 
difficulty of dealing with a continuous spatial variable. Such a problem has been solved by 
introducing a proper discretization 5. As a result, it turns out that in the bulk, 5 = 0(1), 
while, inside b, it is 5 = 0(e). The reason for the difference is that in the bulk, nearby 
configurations change their mutual distances only as a result of local instabilities that are of 
order 1. On the contrary, in b, the difference may also grow due to the propagation of uncon- 
trolled perturbations from the boundaries. Accordingly, the fractal dimension is basically 
proportional to the number of lattice points introduced in the discretization processs: in the 
first and second term of the r.h.s. of Eq. (JSjl . one can recognize the contributions arising 
from the bulk and the boundaries, respectively. The latter one has size 0(l/e 2 ), because the 
length of the boundary is estimated to be 0(1/ s). Since in lattice systems there is a natural 
spacing, an extension of this reasoning to that context would lead to a correction term of 
order 0(l/e), to be confronted with the logarithmic correction predicted by Eqs. 13)5(1 . In 
the last section we will briefly discuss the possible reasons of such a discrepancy. 



IV. LYAPUNOV SPECTRA OF OPEN SYSTEMS 



It is several years that the concept of convective Lyapunov exponent has been successfully 
introduced to describe how perturbations spread and grow. This is done by measuring at 
time t the amplitude 5x l (t) of an initially 5-like perturbation (5x*(0) = 5q) and determining 
its growth rate in a frame moving with velocity v [13I ]. 

A c (v) = km ' Wl , (13) 

t^oo t 

where i = vt. 

It is known that whenever the spatial left-right symmetry is not broken, the maximum 
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value for A c (v) is obtained for v — and it coincides with the standard maximum Lyapunov 
exponent. Upon increasing the velocity, the convective exponent decreases and becomes 
negative for v > v c , to indicate that only perturbations moving with a velocity slower than 
a critical velocity can be sustained. 

One might imagine to generalize this procedure, by looking not just at the amplitude of 
a single perturbation but to the volumes spanned by a finite number of perturbations, very 
much in analogy to what done for computing Lyapunov spectra in closed systems. However, 
if we entirely follow the standard approach, we are bound to conclude that all the convective 
exponents that can be associated to a given velocity coincide with the maximal value. The 
reason for this conclusion is that, on the one hand, the finite number of perturbations that 
are followed in time visit a space of increasing (eventually infinite) dimension. On the other 
hand, the existence of a limit Lyapunov spectrum means, as discussed above, that A; is 
a function of p = l/L alone, where / denotes the l-th exponent, L being the length of 
the system. Since in the above setup, the number of exponents is fixed and equal to L, 
while the system size increases with time, one is basically computing an increasingly thin 
portion of the Lyapunov spectrum and, eventually, all Lyapunov exponents become equal to 
the maximum one. Stated otherwise, by implementing the usual closed systems approach, 
without modifications, we compute the invariant Lyapunov spectrum, A(p), restricted to 
a range < p < p max = L/L e ff, where L e ff is the effective dimensionality of the space 
explored, which increases indefinitely with time, thus implying that p max — > when t — > oo. 

Although this is an inescapable conclusion, we now show that meaningful results can be 
obtained even if perturbations are followed for a finite time. A priori, one might think that 
computing over a finite time implies that the corresponding quantity is ill-defined, because it 
would depend on the choice of coordinates. However, in so far as time is finite but arbitrarily 
large, this objection does not apply. This is for instance the case of the so-called multifractal 
analysis of low-dimensional chaos. In the present context there are two scaling parameters to 
deal with: the length L and the time T. We eventually want to let both diverge to infinity. 
The standard approach adopted in the literature consists in first letting T diverge to infinity 
(this allows determining the Lyapunov spectrum of a finite system) and then taking the 
thermodynamic limit L — > oo. 

The problem of choosing the most appropriate order in the problem at hand is very similar 
to the problem mentioned in the introduction about the order of the two limits L — > oo and 



12 



e — > for a meaningful definition of fractal dimension in open systems. We propose here to 
let T and L diverge simultaneously, with fixed ratio, 

9 = T/L ■ (14) 

i.e., given a subsystem of length L (embedded in a formally infinite chain), we let pertur- 
bations evolve for a time T = gL. Our claim is that the corresponding spectra converge, in 
the limit L — > oo, to a specific shape that depends only on g. 

In order to be more precise, we start defining all the quantities of interest. Consider 
L independent vectors {u n (t)}, (n = 1,...,L). Each u n denotes a perturbation initially 
restricted to a subchain of length L of our (virtually) infinite chain: (that is, Vn = 1, . . . , L, 
it is Un\o) = for i < and i > L, where Un stands for the z-th component of the n- 
perturbation vector). Let us introduce the projection operator, Pl 

r ._ (i£\t) if 0<i<L 
{P L u n (t)}W = I (15) 

I otherwise . 

Let then evolve the L vectors {u n (t)} up to a time T = gL, and determine the volumes 
spanned by the projection of k such vectors (with k — 1, . . . , L), as one normally does in the 
computation of standard Lyapunov exponents. In this way we can compute the Lyapunov 
spectrum over a time T in a spatial window of size L. 

This approach applies to any one-dimensional system, irrespective of the continu- 
ity/discreteness of the space and time variables. However, for the sake of computational 
simplicity, we shall restrict ourselves to consider coupled map lattices. In particular, we will 
mainly refer to the typical coupled map lattice with diffusive coupling, 

x\t + 1) = / [(1 - 2e) x\t) + elx^it) + x i+1 (t)]] ; (16) 

where e is the coupling constant and f(y) is a map of the unit interval onto itself. 
In particular, we start considering Bernoulli maps, where 

f(y) = a(y-[y}) (17) 

and the square brackets denote here the integer part. 

In Fig.^ the spectra corresponding to the same ratio g = T/L = 1, but different lengths, 
have been plotted. One can clearly see a convergence towards an asymptotic limit. The large 
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deviation observed in the bottom part of the solid curve is due to numerical inaccuracies. 
Indeed, these spectra, that are obtained by letting the perturbations evolve in the whole 
space and then projecting them onto the window of interest, require an increasing accuracy 
for increasing elapsed time. For L = 192 even FORTRAN extended accuracy (equivalent 
to approximately 30 digits) is no longer sufficient. Anyway, the inset clearly confirms the 
tendency to converge towards a well defined asymptotic shape, so that one can meaningfully 
introduce the concept of open-system Lyapunov spectrum X(p, g) (OSLS). 




FIG. 1: Open-system Lyapunov spectra of Bernoulli maps with a = 1.2 and g = 1. Dot-dashed, 
dotted, dashed and solid curves correspond to L = 24, 48, 96 and 192, respectively. The last part 
of the solid curve is truncated because it is numerically unreliable. In the inset, an enlargement of 
the region around the maximum is reported. 

The spectra reported in Fig. El correspond instead to different values of g for a fixed 
length L. In the limit g — > 0, the effects of propagation outside the initial window are 
negligible, so that A(p, g = 0) reduces to the standard Lyapunov spectrum. In fact, the 
lowermost curve corresponds to the analytically known expression for the standard Lyapunov 
spectrum. Upon increasing g, the OSLS increases and for g — > oo, we expect it to flatten 
around the maximum Lyapunov exponent. In the inset of Fig. we have suitably rescaled 
the p axis. The rather good data collapse suggests that the propagation of perturbations 

14 




FIG. 2: Open-system Lyapunov spectra of Bernoulli maps with a = 1.2 and L = 96. From bottom 
to top, the curves correspond to g = 0, 1/4, 1/2, 1, and 3/2, respectively. In the inset, the same 
spectra are plotted after rescaling the p axis according to the argument discussed in the text. 

does not modify the spectrum structure, but simply leads to an expansion of the p scale. 
Indeed, from the definition of p = l/L, one can notice that the scale is controlled by the 
window length L. 

However, as anticipated above, in the case of open-system simulations, the space covered 
by each perturbation increases with time, so that it is reasonable that the label I should 
be more properly scaled to some effective length L e // rather than to the initial length 
L. Furthermore, it can be conjectured that the effective length increases linearly in time 
as L e ff = L + 2vT, with some velocity v (the factor 2 is included to account for the 
growth on both sides of the window). Therefore, the "right" scaled variable should be 
p' = l/[L + 2vT] = p/(l + 2gv), and, asymptotically, one expects that X(p,g) converges to 
some A(p'). The good data collapse observed in the inset Fig. El points in this direction, 
although, from the raw data we cannot exclude that the apparent hyperscaling is only an 
approximation. Indeed, the finite size corrections affecting the various spectra are of the 
same order as the deviations among the various curves. 
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In order to study the dependence of the OSLS on g in a more quantitative way, we proceed 
as follows. Given a spectrum X(p,g), we fix a threshold A s and compute the quantity: 

S(X s ,g) = J(X(p,g)-X s )dp, (18) 

where the integral is restricted to the interval of p-values where the integrand is positive. 
For X s = 0, S reduces to the well known Kolmogorov- Sinai entropy. The dependence of 5* on 
g can be observed in Fig. |3] for three different thresholds. The essentially linear behaviour 
confirms that the increase is due to a propagation process, since the effective length increases 
linearly with time, i. e., with g. If only one propagation velocity v s is present in the evolution, 
each curve should increase as 

S(X s ,g) = S(X s ,0)(l + 2v s g) . (19) 

Therefore, the curves obtained for different values of X s should, after rescaling them to the 
same starting point, overlap. That is, if the velocity v s does not depend on the threshold X s , 
then the plot of S'(g) = S(X S , g)/S(X s , 0) should be universal. The inset of Fig. El reveals 
however a weak dependence on X s : the slopes corresponding to the different thresholds are 
0.23, 0.25 and 0.28, indicating that v s ranges in the interval [0.11 — 0.14]. Whether the 
fluctuations are due to finite size corrections or are an indication of a whole spectrum of 
velocities, we cannot say. The increase of the velocity when X s decreases suggests, however, 
that the less unstable directions are characterized by a more efficient propagation. 

Bernoulli maps certainly represent the simplest model for testing the spreading properties 
of perturbations in chaotic systems, but the pecularities of the model (above all, the absence 
of fluctuations for the local multipliers) may invalidate the generality of the conclusions. 
Though this is not true for the scaling properties of the standard Lyapunov spectrum, it is 
nevertheless wise to repeat the analysis above for different models. 

For this reason we now discuss the case of asymmetric tent maps, for which 

V - if < y < b 

b ~ y ~ 

m = < (20) 

l — I if b<y<l. 

While the simulations on this model have been performed for several choices of the 
coupling e, and map parameter b, the results here presented refer mainly to the (e, b) = 
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0.5 1 1.5 g 2 

FIG. 3: The pseudo entropy S(X s ,g) defined in Eq. (|T%|) versus g for A s = —0.7, —0.5 and —0.3 
(from top to bottom), for Bernoulli maps with a = 1.2 and L = 96. In the inset the same data are 
plotted, with the same symbols, after rescaling, giving the quantity S '(g) defined in the text. 

(1/3,3/4) case. Once more, we observe that the OSLS converges to an asymptotic shape 
when simulations are performed for increasing time (and thus window size) at a fixed value 
of g. Some OSLS's are plotted in Fig. where we again see the same flattening tendency 
for the spectrum and a similar collapse, as observed for the Bernoulli maps. In order to test 
quantitatively the (visually suggested) hypothetical hyperscaling, we again investigate the 
behaviour of the pseudo entropy S(X S , g) for various choices of the threshold \ s . The results 
plotted in Fig. confirm what found before, though indicating a more evident increase in the 
slopes (in the rescaled representation) of the fits as A s is lowered, thus reinforcing the idea 
that different velocities are into play 2 . As before, decreasing the threshold X s , from +0.2 to 
0.0 and to —0.2, we observe an increase in the velocity, v s = 0.10, 0.12, 0.14, respectively, 
in agreement with the interpretation, already suggested in the Bernoulli case, of a more 
efficient propagation of disturbances along more stable directions. 

2 The simulations performed with a weaker coupling, e — 1/6, show a qualitatively similar behaviour, 
though, as expected, with a slower propagation, v$ = 0.05 ± 0.01. 
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FIG. 4: Open-system Lyapunov spectra of asymmetric tent maps, for a window of size L = 96, 
with coupling e = 1/3 and b = 3/4. From bottom to top, the curves correspond to g = 0, 1/4, 3/4, 
1, and 3/2, respectively. In the inset, the same spectra are plotted after suitably rescaling the p 
axis. 

Finally, as a prototype of a model with more than one variable per lattice site and also as 
an example of a system with conservation of volumes, we have studied coupled symplectic 
maps, 



where T l (t) = sin [q l (t) — q l 1 (t)]. Consequently, for the evolution of perturbations, we have, 



p\t + 1) = p\t) - K [sing* (t) + e (T l+1 (t) - T(t))] 



q \t + l)=q\t)+p t {t+l) , 



8p\t + 1) 



Sp^t) - K [cosg l (t) - e{D l+1 {t) + D\t))] 5q\t) + 
eK [D l+1 (t)5q t+1 (t) + D i (t)5q i ~ 1 (t)] 



and 



5q\t + l) = 5q l (t) + 5p\t + l) , 



where D l (t) = cos [q l (t) 
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FIG. 5: The pseudo entropy S(X s ,g), defined in Eq. (|18|). for the same tent maps of the previous 
figure, versus g. From top to bottom, it is A s = —0.2, 0.0 and + 0.2. In the inset the ratio 
S' = S(\ s , g)/S(\ s , 0) is plotted, with the same symbols. 

Again, different values of the parameters have been chosen, though we present here just 
the results related to the case K = 5 and e = 1/6. 

In Fig. IH1 the OSLS in the case of a window of size L = 96 (thus, for 2L = 192 Lyapunov 
exponents) are shown, for several values of g. The inset in the same figure points out 
the good scaling behaviour of these spectra. Performing on them the same analysis as 
before, through the computation of S', we find a disturbance propagation velocity v s = 
0.050 ± 0.005 considerably smaller than in the previous cases, though it should be stressed 
that the coupling strength is now smaller by a factor 2. 

The velocities obtained through the analysis performed in this section are basically the 
volume spreading velocities. As mentioned in the beginning of this section, another velocity, 
v c , can be determined from the zero of the spectrum of the convective Lyapunov exponents. 
In order to compare the two velocities, we have computed also v c , using the procedure 
described in Ref . ^| , to which we refer for more details. Here we limit to report the results 
of the application of such a method to the models discussed in this paper. In Table I, the 
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FIG. 6: Open-system Lyapunov spectra of coupled symplectic maps for L = 96 (i.e., N = 192 
exponents). From bottom to top, the curves correspond to g = 0, 1/4, 3/4, 1, and 3/2, respectively. 
In the inset, the same spectra are plotted after suitably rescaling the p axis. 

two sets of velocities can be mutually compared. For a meaningful comparison with the 
velocity obtained from the corresponding convective spectra, we have reported there the 
value of 2v s , since the effective length implicitly determined when estimating the Lyapunov 
spectrum over a time T, is the average length over such an interval and thus it is half of the 
final length (if the growth is uniform, as we are assuming). Therefore each velocity has to 
be doubled, if we want to make a direct comparison with v c . 

We see that in all the cases, the volume spreading velocity is significantly different from 
(and smaller than) the propagation velocity of perturbations, thus indicating that at least 
two mechanisms exist in spatio-temporal chaos that contribute to the propagation of infor- 
mation. Indeed, a moment's reflection reveals a crucial difference between the two velocities. 
On the one hand, we see that if the local multipliers change all by the same constant fac- 
tor, r, the resulting convective Lyapunov spectrum is shifted by lnr and, as a consequence, 
the velocity v c , corresponding to the crossing point with the f-axis, changes. On the other 
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Model 


2v s 




Bernoulli (e = 1/3) 


0.24 ±0.03 


0.48 ±0.01 


Tent (e = 1/3) 


0.24 ±0.04 


0.90 ±0.01 


Tent (e = 1/6) 


0.09 ±0.02 


0.64 ±0.01 


Symplectic (e = 1/6) 


0.09 ±0.01 


0.83 ±0.01 



TABLE I: Synoptic table containing the velocities determined from the open-system Lyapunov 
spectra and from the zero of the convective Lyapunov spectra for the classes of maps discussed in 
the text. 

hand, the velocity v s remains unchanged. In fact, we believe it is not by chance that v s 
turns out to be approximately the same for the same values of the coupling strength (see 
the table): it measures how the spatial coupling forces an ensemble of perturbations to 
cover all neighbouring directions. Thus, in principle, it can be defined even for a stable 
system. In a sentence, we could summarize stating that v c reflects the (chaotic) features of 
the local dynamics, whereas v s exploits the efficiency of the of the coupling in favouring the 
propagation of information. 

V. SCALING BEHAVIOUR OF THE INVARIANT MEASURE AND CONCLUD- 
ING REMARKS 

In this last section we exploit the knowledge of the scaling behaviour of the OSLS to shed 
light on the corrections to the Kaplan- Yorke formula arising from the spatial coupling. We 
start neglecting the integration over the external degrees of freedom (see Eq. ()12j)). In this 
approximation, we are entitled to use equations Q and (fTU|) . with the warning that now the 
Lyapunov exponents do depend on time. Under the assumption that hyperscaling holds, 
Eq. (JTTJ) still applies, with the length L replaced by its effective value over a time T, 

D KY (T,L) =d(l+2vg)L = dL + 2dvT (21) 

In order to independently check the linear dependence of Dry on g, we have studied the 
behaviour of dxy = Dky(T, L)/ L, by integrating the spectra obtained for different values 
of g, in the case of coupled Bernoulli maps. From the above argument, we expect that 
dxY = ^(1 + 2vKYg), where vky is, in principle, still another velocity of propagation, 
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though much of the same nature as v s . 

The data plotted in Fig. [7| indicate indeed a quite clean linear growth. The deviation of 
the value corresponding to g = 1/4 is certainly due to finite-size corrections, since it has been 
obtained for the shortest time compared to the other points, while the datum corresponding 
to g = (d KY {0) = 0.39) is obtained from the standard Lyapunov spectrum (see Ref. 
and is therefore not affected by numerical convergence problems. From the figure, the slope 
of the curve dxyig) is 0.09 and from this (and from cIky(0)) it follows that vky = 0.12, well 
consistent with the value, v s , obtained from the generalized KS-entropy approach. 

We can now use Eq. to transform the dependence of Dky on time T into a dependence 
of D KY on the resolution e. Indeed, inverting Eq. (jUJ), we have 

m hie , „ 

T = — , (22) 

where = X(p = d). 

Inserting then Eq. (|2*2*|) in Eq. ()21j) . we obtain 



which again reveals a logarithmic dependence on the resolution. At variance with the past 
derivations of similar formulas, however, the multiplicative factor in front of the logarithmic 
correction follows now from a well documented discussion of the dynamical evolution. 

However, we should not forget that the above formula represents a lower bound on the 
dimension, Dky, since we have neglected the integral over in Eq. (fT2*|) that accounts 
for the uncertainty on the inner variables induced by the initial lack of knowledge on the 
outer degrees of freedom. An accurate quantification of the corresponding propagation of 
information is far from trivial, since it may be controlled by nonlinear mechanisms as it is 
known to happen in some situations. A typical case, where this is certainly true is that of 
the so-called stable chaos ^(|, i.e. of the irregular behaviour emerging even in the presence 
of a negative Lyapunov spectrum. In that case, no extensive contribution to D exists since 
d = and the previously discussed logarithmic correction is absent too. Nevertheless, finite- 
amplitude perturbations may enter the window of interest carrying relevant information. 

In order to complement the lower bound above, we present here a heuristic argument 
that allows us determining also an upper bound, through an estimate of the effect of the 
disturbances propagating from the boundaries into the region of observation. It should 
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FIG. 7: The scaled dimension cIky versus g, as computed from the Kaplan- Yorke formula, for the 
Bernoulli coupled map lattice. 

be clearly stated that the validity of the following argument relies on the assumption that 
nonlinear mechanisms are negligible. 

From the convective Lyapunov exponents, we know that a perturbation originating at 
the boundary at time 0, is amplified by a factor exp A C (£/T)T, after a time T at a distance 
I. Let us assume that, according to Eq. (}2*2*|) . the time T corresponds to a resolution lne/Arf. 
Then, the maximum distance travelled by the external perturbation while still being larger 
than e can be implicitly determined by the equation 



Let us now v' denote the velocity such that A C (V) = A^ (if the whole convective spectrum 
is larger than A^, v' must be assumed equal to the maximal possible velocity, which, in the 
case of nearest neighbour coupling, is just equal to 1). Eq. (J2H) implies that 

£=^\ne . (25) 
Ad 

By assuming that the distribution of the state variable is independent in all the I sites, one 
finds that f£ (where / denotes the number of variables per lattice site), represents an upper 




(24) 
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bound to the correction due to the external degrees of freedom. Accordingly, this correction 
has again the same structure as the previous one; altogether, one finds that 

D c {s,L)<dL + 2i ~ dV ^ fV ' ) Ins. (26) 

Thus, we see that the dependence of the effective dimension on e remains weaker than the 

n 

l/e behaviour deduced indirectly from the treatment in Ref. |4|. We interpret this as the 
indication that the rigorous bound determined in Ref. j^J can be improved. Neverthless, a 
more detailed understanding of the propagation of perturbations is also required, in order 
to shed further light on the exact structure of the logarithmic corrections. We are currently 
exploring the possibility to directly quantify the amplitude of volume-perturbations in a 
simplified class of systems where most of the calculations can be carried on analytically jl7|. 

The above analysis has been perforemd under the implicit assumption that |lne| < 
aL (where a is a suitable dimensional factor). This limitation is equivalent to the one 
encountered in the paper by Collet & Eckmann [4], where it is stated that the bound (jBJ) 
applies only in the region where e > a/L. The reasons for these limitations are due to 
nonlinear effects. In the previous section we have seen that in the linear approximation 
any perturbation, initially restricted to Wl, sooner or later diverges along all the L existing 
directions (inside the window). However, nonlinear corrections may become important much 
before the amplitude of the perturbation becomes of order 1 inside Wl- Indeed, we have 
seen that perturbations grow outside Wl, too; in particular, those that initially decay inside 
Wl, grow outside the window. As soon as their amplitude becomes 0(1) in the outer part, 
the effect of nonlinearities cannot be any longer neglected even inside Wl, though their 
amplitude is therein very small. Since we have seen that high resolutions correspond to 
relatively long time, we cannot expect our approach to hold for those tiny scales where 
| In e | /L is no longer small. 

It is now worth commenting about the relationship between the open-system Lyapunov 
exponents introduced in this paper and the previously devised chronotopic Lyapunov ap- 



proach 



18, 



19l | . There, it was conjectured that all linear stability properties of Id spatially 



extended dynamical system can be obtained from the entropy potential: a function of the 
spatial and temporal growth rate of generic perturbations. Although we have not been able 
to find the specific link with the OSLS, there is no reason to think that the information con- 
tained in that class of Lyapunov spectra is not contained in the entropy potential. Finding 
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the relationship between the two approaches would be very interesting not only from a con- 
ceptual point of view, but also because it would allow for a much easier computation of the 
OSLS. We must, indeed, recall that an accurate determination of Lyapunov spectra in open 
systems is hindered by the high accuracy that it requires: this limitation is certainly crucial 
in the context of continuous space-time systems such as, e.g., the complex Ginzburg-Landau 
and the Kuramoto-Sivahsinsky equations. 

Finally, we expect the open system approach to be of some relevance also in connec- 
tion with the nonequilibrium dynamics of Hamiltonian systems. For instance, in Ref. j^j], 
Gallavotti was interested in the expansion rate of local (in space) volumes. This is nothing 
but the function S, defined by Eq. (JTHJ), for \ s corresponding to the minimum Lyapunov 
exponent. 
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